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1. Introduction 

The strong coupling constant a s = g 2 / (4n) is one of the Nf+l input parameters of QCD and 
as such one of the fundamental constants of nature. Its actual value depends on both the renormal- 
isation scheme (including the number of active flavours) and the scale. Given a renormalisation 
scheme S the dependence on the scale pL is controlled by the renormalisation group (RG) through 

2 d a? (a 2 ) oS , SN of-j-o os f a f 
d\l l % gj V 71 

where /3 s is the beta function in that scheme. Solving Eq. ( |Ll| ) yields an exact relation between 
the scale-dependent coupling of (/i) and the RG-invariant, scale-independent but renormalisation- 
scheme-dependent parameter A s , defined via ln/x 2 /Ag = / daf /j8 s (a). Once the Lambda param- 
eter is known for one scheme, a one-loop calculation suffices to determine it in any other scheme. 

While recent precision determinations of a s exist, based on either perturbative analyses of 
short-distance-sensitive lattice observables or sum rule analyses of hadronic z decay data (for de- 
tailed discussions and other relevant references, see [[I], ||] and, e.g., ||3|], respectively), residual 
uncertainties mean that additional independent high-precision determinations remain of interest. 
Some of us have recently introduced the minimal momentum subtraction (MiniMOM, or MM) 
scheme for QCD in covariant gauges [|J]. An important advantage of this scheme is that it allows 
the strong coupling to be fixed solely through a determination of the gluon and ghost propagators. 
In Landau gauge this scheme has been implicit in the early studies of these propagators []|]. 

The MM scheme is defined by combining MOM scheme propagator renormalisation with the 
supplementary condition Z\ = Z MS for the ghost-gluon vertex renormalisation constant With 
Z and G the respective gluon and ghost dressing functions, the MM coupling is then defined as ^ 

aT(P 2 ) = ^Z(p 2 )G 2 (p 2 ). (1.2) 

The relation between a s MM and a ? MS is known to four loops PJ]. Here we use a s MM to determine 
A-j^ (in units of ro) for Nf = 0,2 from continuum extrapolations of the product of the bare lattice 
Landau gauge propagators, as first proposed in [Q]. 1 The absence of vertex measurements in the 
method allows for a significantly improved accuracy in the lattice estimate for a 4 MM . 




2. Numerical setup 

The results below were obtained on both Nf = and 2 SU (3) gauge field configurations. The 
quenched configurations were thermalised using the standard Wilson gauge action, with j8 ranging 
from 6.0 to 8.5, applying standard update cycles of heatbath and micro-canonical over-relaxation 
steps. The unquenched gauge field configurations were provided by the QCDSF collaboration, 
who used the same gauge action supplemented by Nf = 2 clover-improved Wilson fermions at 
various values of the hopping-parameter K (see Tab. |l| for further details). All gauge configurations 

'Our 4-loop expansion for j3 MM (o!) was not worked out until January 2008. While the 3-loop version of this 
expansion in fact differs somewhat from the 3-loop MOM/; scheme result used in uh, the difference is small. 
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Table 1: Parameters of j 
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were fixed to lattice Landau gauge using an iterative gauge-fixing algorithm. To guarantee high- 
precision the local violation of transversality was not allowed to exceed e < 10~ 10 where, as usual, 
£ = max* <KeTr [(V A1 A. v )(V^A XAi ) t ] and A xfl = ^(U x ^ - tfJ M )ltraceiess ■ 

Gluon and ghost propagators were measured on these gauge-fixed sets employing standard 
techniques and an acceleration for the Faddeev-Popov-operator inversion (see ^] for details). 
As verified numerically in 0, for the range of momenta studied here the Gribov ambiguity is 
irrelevant. Using well established values for ro/a [|8]-[l0|] to bring the raw data on a^ M for different 
j3 (and fc) onto the common scale r\p 2 (see Tab. [l]for the ro/a values), and with g 2 (a) = 6//$(a) 
the bare coupling at the lattice cutoff scale a -1 , a s MM was then determined from the averaged data 
for the bare lattice gluon and ghost propagator dressing functions, Zl and Gl, via 

a s MM (p 2 ) = « L MM (p 2 ) + 0(a 2 ) with ar(p 2 )^ 8 -^Z L (p 2 ,a 2 )G 2 L (p 2 ,a 2 ). (2.1) 

To have the lattice tree-level structure correct the dressing functions were extracted using aq^ ) = 
2sin(a/7^/2), but a s MM is considered versus p 2 with ap^ = Ink^/L^ and G (—L^/2,L^/2\. 

3. Lattice data of the MiniMOM coupling 

Our lattice data for a^ M (p 2 ) is shown in Fig. |], together with the expected 4-loop continuum 
running obtained using, to be specific, the ALPHA collaboration MS values, *"oA^ = 0.60(5), and 

r A^ = 0.62(4) (4) [JllJ, translated to the MiniMOM scheme (see Table 4 of Ref. g] for the 
relevant values of A^/Amm)- m what follows, for the sake of illustration, we display momenta 
in physical units (rather than as r\p 2 ), using ro = 0.467 fm = 2.367 GeV -1 . Since only ratios of 

(Nf) 

momenta enter the calculation, this choice does not affect the final result for rnA-^L . 

MS 

From Fig. [j], one sees that, for both Nf = 0, 2, scaling violations, finite volume effects and 
hypercubic lattice artefacts are nearly negligible, even though, for Nf = 0, the lattice spacing varies 
over an order of magnitude. Small systematic deviations from continuum 4-loop running, however, 
become visible at higher resolution. Such deviations are negligible for purely diagonal lattice 
momenta satisfying 3 < a 2 p 2 < 6, but grow, oc \/{ap) 2 (c< a 2 p 2 ) to leading order, for momenta 
below (above) this interval. Fig. [j] shows data for diagonal momenta with 1 < a 2 p 2 < 10. 

Deviations from continuum 4-loop running are not unexpected. At small momenta, they result 
from a mixture of (a) finite volume effects, (b) the onset of nonperturbative effects (condensates 
etc.) and (c) truncation errors in the perturbative expansion of the coupling, while at large momenta 
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Figure 1: Lattice data for a^ M (p 2 ) for Nf = (top) and Nf = 2 (bottom) compared to the expected 4-loop 
running of af AM (p 2 ) (lines) in the continuum limit Note that these are not fits: The data is the raw 
lattice data for purely diagonal momenta ap^= Ink^/L^ with 1 < a 2 p 2 < 10 brought onto a common scale 
employing established values for ro/a for these sets [^|-|l0||. For illustration purposes, the 4-loop running 
has been fixed using r Q A^ = 0.60(5) and r A^ = 0.62(4) (4) of the ALPHA collaboration [JTT|, ||], and 
the overall momentum scale set using ro = 0.467 fm = 2.367GeV~ 1 . 
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they are due to (d) scaling violations proportional to a 2 p 2 and (e) hypercubic lattice artefacts pro- 
portional to higher-order invariants a n p^ = £^ a n p n ^ (n = 4, 6, 8) of the isometry group H(4). The 
latter, in particular, would become pronounced for larger momenta without suitable corrections, 
which can be performed either using the H4 method or by imposing so-called cylinder cuts on the 
data (for the former approach, see, e.g., [|i"3|], for the latter Jl4|]). The cylinder cut approach, though 
less sophisticated, is nonetheless effective and robust. We chose a combination of the two methods 
to keep the statistical noise to a minimum. To be specific, we consider data only for purely diagonal 
lattice momenta, for which hypercubic lattice artefacts are known to be smallest, and correct for 
the remaining (rather small) artefacts through a fit of this data to a hypercubic Taylor expansion of 
the (lattice) M M coupling which, to leading order, has the form 

arip 2 ) = oc^ M (p 2 ){l + c l .a 2 p 2 + c 2 .aY + ...) (3.1) 

where the c,- are constants (see also [jlj, [TJ], [l6|]). 
For the classic H4 method one would have to fit 
the C;'s from extrapolations of the a^ M (p 2 ) data 
at different lattice momenta belonging to differ- 
ent H(A) orbits but the same a 2 /? 2 . Not only are 
these extrapolations susceptible to statistical arte- 
facts for insufficiently large lattice sizes, but the 
Faddeev-Popov operator has to be inverted us- 
ing point sources to get data for all momenta, 
introducing larger statistical fluctuations into the 
ghost propagator (and hence the coupling) at large 
momenta. Our approach allows us to use instead 
plane-wave sources for our inversions (see ^ for 
details) and fully exploits the translation invari- 
ance of the lattice, thus drastically reducing sta- 
tistical noise in the coupling. The results below 
bear out the reasonableness of the approach, providing an excellent description of the data for 
3 < a 2 p 2 < 30 (or 3 < a 2 /? 2 < 12 for C2 = 0) and giving stable results for all fit parameters. 

The deviations at small momenta, illustrated in Fig. § are more severe and, in our opinion, not 
yet fully under control. They start to become visible for a 2 p 2 < 1 and appear to be a mixture of finite 
volume and nonperturbative effects (the latter expected to set in at smaller j8). The effect is such 
that data points at fixed physical momenta decrease as the physical volume increases. Currently, 
additional simulations at different j3 but fixed physical volume, are being performed to help bring 
these low-momentum artefacts under better control. For now, we exclude data with a 2 p 2 < 3 to 
stay well clear of the region where such effects become evident. 




10 100 1000 

p 2 [GeV 2 ] 



Figure 2: Nf = data for a™ M (p 2 ) as shown in 
Fig. 0but here for < a 2 p 2 < 23 to illustrate finite 
volume effects at small momenta which are more 
pronounced here than the lattice artefacts at large. 



4. Fitting the data 

Our fitting procedure works as follows. Each of the data sets is fitted separately to the 
Ansatz (3.1), where the 4-loop perturbative running form is used for a s . (j> 2 ) and the remain- 
ing terms correct for the leading lattice artefacts at larger a 2 p 2 . All fits are performed using the 
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Figure 3: Values for Q^ f s obtained from fits to the data for different afro using Ansatz (3.1 ); left for Nf = 
and right for Nf = 2 and different bare quark masses amo. See text and Tab. [l]for further details. Symbols 
refer to the same /3, N 4 (and k) as in Fig. [j]. The dotted line and very left 'star' (left panel) are a fit of the 
plateau of the three points (filled symbols), calculated in equal physical volumes V[ « (1.2fm) 4 . Note that 
there we have also included two very recent points (brown crossed circles) from simulations at /3 = 6.7 and 
6.92 and equal physical volume, V2 ~ (1.6fm) 4 . 



fixed fitting window 3 < a 2 p 2 < 12. (The fits have been checked to be quite robust to small changes 
to the lower and upper bounds of this window.) We also find that C2 can be set to with little effect 
on c\ . The 4-loop running of a 4 MM (p 2 ) is generated from that of a s MS (p 2 ) using the 4-loop relation 
between the couplings given in Eqs. (14-15) of Ref. [Q]. (Further, more specific details will be 
provided in an upcoming publication.) 

With C2 set to zero, our fit parameters reduce to c\ , the constant of the leading lattice correction 
at larger momenta, and 0^ f s , the MS coupling at an arbitrary reference scale p 2 ei (to be specific, 
we take p 2 ei = 70GeV and ro = 0.467 fm). A^g could, of course, be used in place of a r ^ f s as a 
fit parameter; we expect our fits to be more stable with the latter choice. A^g in any case follows 
from using the standard relation [17] 



jj. 2 (da 1 f 1 



A 2 J P(a) ft 



^+b l lna+(b 2 -b 2 l )a+ fy-£i& 2 + y ) a 2 



+ C 



where a{p.) = a s {jx)/n, j6 ( ^ s , . . . ,j8 3 MS are the MS scheme j3 -function coefficients, b t = fr/po and 
C = (Vj3o)lnj6 . 

Fitted values for aM s and ci as a function of a/ro (and, for Nf = 2, also for different bare 
quark masses) are shown in Figs. || and |j. Note that if the data showed perfect scaling, all fits 
would give the same value (within errors) for a r ^ f s , independent of the lattice spacing, and the 
fitted 0(a 2 ) corrections (ci) would turn out to be zero. As expected, c\, though small, is not zero 
(see Fig. |]). The long plateau for c\ as a function of a/ro, however, suggests that our Ansatz for 
describing the leading lattice corrections is a reliable one. This conclusion is also supported by the 
results for Of r ^ f s . For Nf = 0, for example, a r ^ f s starts to plateau around a/ro = 0.1, i.e., for /3 > 6.4 
(see Fig. ^). Note that the fit quality is significantly degraded if no correction term is included. 

For Nf = 2 we observe only small deviations in a^ f s on changing the bare quark mass (see 
Fig. ||, right panel). Changing a/ro leads to more significant shifts. This is almost certainly due 
to the relatively large a/ro employed, and the fact that the lattice data available for a A MM come 
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Figure 4: Coefficient c\ of the leading lattice correction (see Eq. (p.l[)) for different a/ro, obtained from the 
same fits as a r ^ f s shown in Fig. || Symbols are the same as in Figs, [j] and ^. Dotted lines mark the average 
c\ including data left for a/ro < 0.15 and right for a/ro < 0.2. 

from momenta where deviations from the 4-loop running of q s M have already set in (see Fig. [I], 
bottom panel). The overall picture, however, appears similar to that for Nf = 0, and the fitted a r ^ f s 



for the j6 = 5.40 configurations are already close to what is expected, for example, from [12]. In 
fact, r A^ = 0.62(4) (4) corresponds to affi = 0.150(3) in the right panel of Fig. |. 

A continuum extrapolation in our approach corresponds to fitting data plateaus visible at small 
enough a/ro. While it is far from clear that a plateau has been reached for the Nf = 2 data in Fig. || 
the Nf = a r ^ f s results level off quite nicely at smaller a/ro. Fitting the three points having the same 

„MS 



physical volume, i.e., those with (fi,N) = {(6.6,32), (6.9,48), (7.2,64)}, yields oQ b = 0.131(1) 

MS 



This corresponds to roA^ = 0.62(1) which agrees well with values from the literature. Assuming 



that the relative decrease of a r ^ f s from a/ro =0.15 to ~ 0. 1 will be the same for Nf = and 2 (about 
1.5%), one would expect the Nf = 2 data to level off at smaller a/ro around a r ^ f s = 0.148(2). 
This would correspond to ^oA^ = 0.59(3), also in good agreement with existing values in the 
literature. Additional data for Nf = 2 at smaller a / ro would allow us to make further progress, but 
adequate gauge configurations are unfortunately not yet available. For now we take the average 

(2) 

r A¥L = 0.60(3) (2) with the second error due to the uncertainty in the continuum extrapolation. 
5. Conclusions 

In this paper, we have taken advantage of the recently introduced MM scheme for QCD in 
covariant gauges [Q] to perform a determination of the QCD Lambda parameter for Nf = 0,2. The 
scheme allows the strong coupling constant, and hence A^g-, to be determined from measurements 
of ghost and gluon two-point functions on the lattice. The restriction to measured two-point func- 
tions, and the fact that the relation between the MM and MS couplings is known to 4-loop order, 
allows for a high precision determination, with reliable error estimates. 

Our results to date are restricted to Nf = 0,2, and must be extended to 7Y/ = 2 + 1 in order to 
reach the desired goal of estimating a s MS (Mz). Fortunately, over the last few years, the number of 
available Nf = 2 + 1 gauge configurations has increased significantly, and continues to increase. 
The Nf = 0,2 results presented above, which yield the (still preliminary) results 

r A^ = 0.62(1) and r A^ = 0.60(3) (2), (5.1) 

in agreement with other studies (e.g., [jl^- 12, I(J), demonstrate the reliability and accuracy of our 
method, and thus pave the way for future Nf = 2 + 1 analyses. The analysis also provides valuable 
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information on how to bring lattice artefacts under control. 

A positive feature of the current study is that lattice artefacts are found to be almost neg- 
ligible if one restricts the analysis to strictly diagonal lattice momenta satisfying 3 < a 2 p 2 < 6. 
For larger momenta, hypercubic lattice artefacts become visible. For a 2 p 2 < 12 they grow like 
ci a™ M (p 2 )a 2 p 2 , with ci = -0.00016(3) for N f = 0, and c\ = 0.00023(4) for N f = 2, and thus 
can be corrected for quite precisely. It is anticipated that the precision could be improved further if 
data from lattice perturbation theory was available (see, e.g, [ |l8| ] for steps in this direction). 

For smaller momenta, finite volume effects are present, in particular for a 2 p 2 < 1. These 
effects are not yet fully under control and new calculations at different j3 but fixed volumes are 
under way to rectify this situation. It is for this reason that we have not yet investigated condensate 
effects, which are expected to be relevant for lower momenta (see, for example, JT^]). 

Note that the M M coupling could also be employed to determine the lattice spacing depen- 
dence a(/$) via aAjj^(P). With this information, and the ratios of the A^- for Nf = 0,2 and 3, a 
well-established experimental value for q $ MS (Mz) could be used to fix the physical scale of a for 
different Nf (rather than using tq). This again would require a good understanding of all lattice 
artefacts since q s MM data for different lattice spacings would have to be brought to a common scale 
via a matching procedure, starting in the perturbative region. 
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